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We explore the potential energy surfaces for NH molecules interacting with alkali-metal and 
alkaline-earth atoms using highly correlated ab initio electronic structure calculations. The surfaces 
for interaction with alkali-metal atoms have deep wells dominated by covalent forces. The result- 
ing strong anisotropics will produce strongly inelastic collisions. The surfaces for interaction with 
alkaline-earth atoms have shallower wells that are dominated by induction and dispersion forces. 
For Be and Mg the anisotropy is small compared to the rotational constant of NH, so that collisions 
will be relatively weakly inelastic. Be and Mg are thus promising coolants for sympathetic cooling 
of NH to the ultracold regime. 
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I. INTRODUCTION 

In recent years there has been growing interest in 
the production and properties of cold molecules. Pos- 
sible applications, such as controlled ultracold chemistry 
[l|, quantum information and computing 0, and high- 
precision measurements of the time-dependence of fun- 
damental 'constants' [3, 0, 0], make cold molecules ex- 
tremely interesting across many different fields of physics. 

Two main approaches to the production of cold 
molecules can be distinguished. One approach is based 
on the coherent formation of ultracold molecules such 
as Cs2 or RbCs in trapped ultracold atomic gases [1]. 
The molecules may be formed cither by photoassociation 
Q or by Feshbach resonance tuning Q. They inherit 
the /xK-nK temperatures of the parent ultracold atomic 
cloud and usually need very little further cooling. Efforts 
in this area have led to the Bose-Einstein condensation 
of Feshbach molecules [1, [13, [lH and to the transfer of 
Feshbach molecules to low- lying states [l2|, [ij, . There 
have also been considerable successes in dire ct p hotoas- 
sociation to produce low- lying states [l^, [l^, llTl . [Tsj . 

In the other ap pro ach, represented for example by 
Stark deceleration [1^, [2^ or helium buffer-gas cooling 
[2^ . preexisting molecules are decelerated either by ex- 
ternal fields or by collisions with other particles and 
trapped in electrostatic or magnetic traps. The tem- 
perature of the resulting molecular cloud is usually in 
the K-mK region, and therefore new ways for cooling the 
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molecules further are being sought. A promising route 
to cooling decelerated molecules down to the /iK region 
is offered by sympathetic cooling. 

Sympathetic cooling, in which one species is cooled by 
thermal contact with another much colder species, was 
originally developed as a cooling technique for trapped 
ions [1^. Diatomic [1^ and polyatomic [2^ molecular 
ions have been cooled to sub-Kelvin temperatures by 
thermal contact with cold ions, and sympathetic cool- 
ing is expected to be capable of cooling ions of very high 
mass, including those of biological relevance. Sympa- 
thetic cooling has also been successful in producing ul- 
tracold neutral atoms of species that arc not themselves 
suitable for evaporative cooling; for example it was used 
to create the first Bose-Einstein condensates of potassium 
[13. 

Sympathetic cooling is effective only if the rate of elas- 
tic collisions is very large compared to the rate of inelas- 
tic collisions. Elastic collisions exchange kinetic energy 
between molecules and allow thermalization. However, 
inelastic collisions in which internal energy is converted 
into relative kinetic energy cause trap loss (if the energy 
released is greater than the trap depth) or heating (if the 
energy released is less than the trap depth). Magnetic 
and electrostatic traps always trap molecules in low-field- 
seeking states, which are not in their lowest state in the 
applied field. Since typical traps have depths in the /iK 
or mK range, most inelastic collisions cause trap loss. A 
commonly stated rule of thumb is that the ratio of elas- 
tic to inelastic collision rates must be at least 100 for 
effective sympathetic cooling. 

The obvious coolants for the sympathetic cooling of 
cold molecules are alkali-metal atoms, which can be 
cooled to ultra-low temperatures on demand. Sympa- 
thetic cooling of photoassociated alkali-metal dimers in 
triplet states by alkali-metal atoms was discussed by 
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Cvitas et al. [26|, |27[, who concluded that the dimcrs 
would need to be in their ground rovibrational state be- 
cause of the unfavorable ratio of inelastic to elastic cross 
sections in the sub-mK temperature region. 

Sympathetic cooling of decelerated molecules by alkali- 
metal atoms was first considered by Soldan and Hutson 
[2^ , who studied the interactions of rubidium atoms with 
NH molecules. They showed that the '^A" and '^A" states 
of RbNH (bound by covalent and dispersion forces) are 
crossed by much deeper ^A' and ^A" ion-pair states in the 
energetically allowed region at linear geometries. They 
concluded that the ion-pair states are likely to have im- 
portant consequences for the physics of sympathetic cool- 
ing of molecules such as CH, NH and OH, because they 
provide additional mechanisms for inelastic collisions and 
three-body recombination. 

Lara et al. [2^, [s^l subsequently focused on the inter- 
action of OH molecules with ultracold Rb atoms. They 
developed full sets of coupled potential energy surfaces 
and carried out quantum collision calculations including 
spin-orbit and hyperfine coupling. Once again they found 
a deep ion-pair state {^A' for RbOH) that crossed the co- 
valent states at energetically accessible geometries. How- 
ever, even when the ion-pair state was excluded from the 
calculation, the anisotropy of the potential for the cova- 
lent states was enough to cause strong inelastic collisions 
that would prevent sympathetic cooling except for atoms 
and molecules in their absolute ground states. Lara et 
al. concluded in general that (i) light atomic partners 
are desirable as coolants, because the resulting high cen- 
trifugal barriers would suppress many inelastic channels; 
(ii) weak coupling of the electron to the axis, which oc- 
curs in Hund's case b molecules such as NH or CaH, 
would be beneficial; (iii) the anisotropy of the atom- 
molecule surface should be comparable to or smaller than 
the rotational constant of the molecule; and (iv) closed- 
shell coolants, such as alkaline-earth atoms, may produce 
more isotropic potential energy surfaces than open-shell 
coolants, such as alkali-metal atoms. Later a detailed ab 
initio study of 2D adiabatic potential energy surfaces for 
NH interacting with Rb and Cs atoms was reported by 
Tacconi et al. [31]; followed by studies of the quantum 
dynamics of ultra- low-energy collision processes (32l . [33| . 

Another set of potential coolants for sympathetic cool- 
ing are alkaline-earth atoms. Calcium [34| and strontium 
[35| atoms can be cooled and trapped at temperatures of 
the order of /iK. Mehlstaubler et al. [s^ have recently 
succeeded in cooling magnesium atoms to sub-Doppler 
temperatures of 500 fxK. To our knowledge no attempt 
has yet been made to cool Be atoms. 

Very recently, Zuchowski and Hutson [sj surveyed 
interactions of NH3 molecules with alkali-metal and 
alkaline-earth atoms. All the systems exhibited deep po- 
tential wells (890 to 5100 cm~^) when the atom was on 
the N side of the molecule, and shallow potential wells 
(100 to 130 cm~^) when the atom was on the H side, re- 
sulting in very strong anisotropy of the surfaces. This 
will produce strong inelasticity in the molecular rota- 



tional and inversion degrees of freedom and sympathetic 
cooling is unlikely to be successful for molecules in low- 
field-seeking states. 

In this paper we survey the possibilities for sympa- 
thetic cooling of NH molecules, which have very re- 
cently been cooled and magnetically trapped at 0.7 K 
in their ground X'^T,~ state [sl] by buffer-gas cooling. 
NH molecules in their metastable a^A state have also 
been Stark-decelerated [3§| and electrostatically trapped 
at ternperatures of 60-100 mK, and there is a pro- 
posal [401 to transfer a^A molecules to the X^T,~ state. 
In the present paper we investigate the interactions of 
NH(X'^E~) molecules not only with all the relevant 
alkali-metal (Alk) atoms, but also for the first time with 
the alkaline-earth (Ae) atoms. We characterize the po- 
tential energy surfaces of the covalent and dispersion- 
bound states of the AlkNH and AeNH systems and lo- 
cate their conical intersections with the ion-pair states. 
We show that Be and Mg atoms are promising candidates 
for sympathetic cooling of NH(Ar'^I]^) molecules. 



II. METHODS 

To facilitate future quantum dynamics calculations, all 
results are reported in Jacobi coordinates. The inter- 
molecular distance R is the distance between the alkali- 
metal or alkaline-earth atom and the center of mass of 
the NH molecule. The angle 9 is measured at the center 
of mass and is zero for linear atom-HN geometries. In 
all our calculations the NH bond length r is fixed at the 
experimentally determined equilibrium value for the free 
monomer, 1.0367 A [4l[. 

Supermolecular coupled-cluster calculations were car- 
ried out using a single-reference restricted open-shell vari- 
ant [l^l of the coupled cluster method [i^ with single, 
double and non-iterative triple excitations [RCCSD(T)]. 
All electrons from the "outer-core" orbitals (Is, 2s2p, 
3s3p, 4s4p, and 5s5p for Li and Be, Na and Mg, K and 
Ca, Rb and Sr, and Cs, respectively) were included in 
the RCSSD(T) calculations. All the ab initio calcula- 



tions were performed using the MOLPRO package [44 1. 

To describe the interaction between the outer-core 
and valence electrons, and to reduce basis-set superposi- 
tion errors, rather large basis sets are needed. We use 
the correlation-consistent polarized valence quintuple- 
C (cc-pV5Z) basis sets of Dunning [i^ for hydrogen 
(without the g functions) and for nitrogen (without 
the h functions). Both these basis sets were aug- 
mented in an even-tempered manner and used in un- 
contracted form. For lithium, beryllium, sodium, mag- 
nesium, potassium, and calcium atoms, we use the 
correlation-consistent polarized core-valence quintuple-^ 
cc-pCV5Z basis sets of Iron et al. [i^, again without 
the h functions. The Li, Be, Na, and Mg basis sets 
were used in fully uncontracted form and those for K 
and Ca were used partially contracted. These basis sets 
were also augmented by additional even-tempered dif- 



3 



fuse functions. The resulting aug-cc-pCV5Z basis sets 
consist of (15s,10p,8d,6/,4(7), (20s,13p,9(i,7/,5(7), and 
[13s,12p,9(i,7/,5g] functions for Li, Na, and K, respec- 
tively, and of (15s,9p,8d,6/,45), {21s,15p,9d,7f ,5g), and 
[13s,12p,9(i,7/,5(7] functions for Be, Mg, and Ca, re- 
spectively. For rubidium, strontium, and cesium, we 
use the small-core scalar relativistic effective core po- 
tentials ECP28MDF and ECP46MDF [13, Sg, together 
with the corresponding valence basis sets. These ba- 
sis sets were augmented in the even-tempered manner 
and used in uncontracted form. The resulting basis sets 
consisted of (14s,llp,6(i,4/,25), {15s,12p,7d,5f ,2g), and 
(13s,12p,6(i,4/,3(7) primitive Gaussian functions for Rb, 
Sr, and Cs, respectively. 

All interaction energies are calculated with respect to 
the separated-atom- molecule limit, with both the atom 
and the molecule in their ground states. The full coun- 
terpoise correction of Boys and Bernardi [1^ is used to 
compensate for basis set superposition errors. Optimiza- 
tions of the counterpoise-corrected dimer interaction en- 
ergies are performed making use of a general optimization 
algorithm implemented in MOLPRO. 

In order to describe the dispersion-bound state of 
MgNH for all geometries, we use a version of symmetry- 
adapted perturbation theory (SAPT) based on a density- 
functional theory (DFT) description of the monomers. In 
the SAPT(DFT) method ^ (which we use here only for 
MgNH), the interaction energy is obtained as a sum of 
contributions. 
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where E^l^^^ is the electrostatic energy, iJgxch first- 
order exchange energy, E^-' and E^^^^ are the second- 
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are their exchange counterparts. The first- 



and E'^ 

order terms are calculated using Kohn-Sham orbitals, 
while the dispersion, induction and exchange-induction 
terms are evaluated using coupled Kohn-Sham density 
susceptibilities. The exchange-dispersion term is esti- 
mated as described in Ref. |5C| . However, in Ref. (s^l the 
second-order exchange corrections are given in the so- 
called approximation, which neglects terms of third 
and higher powers in the overlap matrix S. Since the 
overlap between Mg and NH is large, we scale the second- 
order exchange corrections by £^oxcii('5'^)/^ciLj This 
procedure was introduced by Patkowski et al. [51| to im- 
prove the performance of SAPT for systems with very 
diffuse monomer densities. 

Our SAPT calculations use the PBEO density func- 
tional [s^ . with augmented correlation-consistent po- 
larized valence quadruple-zeta (aug-cc-pVQZ) basis sets 
[45I . l46j supplemented with bond functions 0.9,0.3,0.1 s 
and p, 0.6,0.2 d and / placed at the midpoint between the 
Mg atom and the center of mass of NH. The Tozer-Handy 



asymptotic correction [53] of the exchange-correlation 
potential is used, with splicing parameters 3.5 and 4.7 A. 



III. RESULTS AND DISCUSSION 
A. Alkali-metal atom + NH interactions 

In discussing the electronic structure of NH interact- 
ing with alkali metals (Alk), it is convenient to begin 
with linear arrangements. At linear geometries (point 
group Coov), there are two covalent states and 
which correlate with the Alk(^S') -I- NH(^S~) dissocia- 
tion limit. These are crossed by an ion-pair ^H state, 
which in a diabatic representation correlates with the 
Alk+(iS') -I- NH"(2H) dissociation hmit. (In an adia- 
batic representation this state changes configuration at a 
long-range avoided crossing with a higher ^H state, and 
in the new configuration it correlates with Alk(^P) -|- 
NH('^S^); the ion-pair configuration is then carried up 
by a cascade of similar avoided crossings until it reaches 
the Alk+(iS') + NH-(2n) dissociation hmit). At non- 
linear geometries (point group C^), the and 
states become ^A" and "^A" states, while the ^H state 
is subject to the Renner- Teller effect and splits into two 
states with the electron hole either in the triatomic plane 
{^A') or perpendicular to it i^A"). In cuts through the 
potential at fixed N-H distance, the covalent ^A" and ion- 
pair '^A" states avoided-cross at nonlinear geometries but 
form a conical intersection at linear geometries. In the 
full three-dimensional picture they form a seam of conical 
intersections parameterized by the N-H distance. 

We consider first the quartet states. The potential 
curves for the states at linear Alk-NH geometries 
are shown in Fig. [1] Results from geometry optimization 
are given in Table H] for both linear configurations. 
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FIG. 1: One-dimensional cuts through the lowest E~ poten- 
tial energy surfaces of AlkNH systems at the linear Alk-NH 
arrangement. 



From the viewpoint of sympathetic cooling, the 
most important feature of the potential surface is its 
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TABLE I: Lowest ''E and ^11 states of linear AlkNH: minima (-Rmin, Knin) and crossing points (_Rx, V^) at difTerent arrange- 
ments (Alk-NH, Alk-HN). Energies are given in cm~^ and distances in A. 
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mm 


Li 


2.176 


-1799.1 


3.09 


-600 


1.78 


-21073 


4.658 


-115.3 


Na 


2.737 


-651.3 


3.21 


-470 


2.13 


-13287 


4.872 


-98.9 


K 


3.073 


-784.7 


3.84 


-432 


2.42 


-14235 


5.386 


-91.1 


Rb 


3.254 


-709.3 


3.99 


-412 


2.53 


-13918 


5.521 


-87.3 


Cs 


3.435 


-737.9 


4.31 


-381 


2.65 


-15116 


5.761 


-85.0 



anisotropy. For sympathetic cooling to be successful for 
low-field-seeking molecular states, we need a potential 
surface where the anisotropy is small compared to the 
rotational constant of the monomer (i?e = 16 cm~^ for 
NH). The global minimum for quartet states is at lin- 
ear Alk-NH geometries for all systems. For LiNH, the 
well depth at the global minimum is about 1800 cm~^, 
while at the secondary minimum (Li-HN) it is only 115 
cm~^. The absolute well depths for the other AlkNH 
systems are a factor of 2 to 3 smaller than for Li-NH, 
while the depths of the secondary wells are comparable 
for all the alkali-metal atoms. Nevertheless, for all the 
AlkNH systems the anisotropy is very large compared to 
the rotational constant of NH. 

The large anisotropy of the AlkNH quartet surfaces re- 
sults from strong sp mixing of the alkali-metal orbitals in 
the Alk-NH arrangement. Such mixing is much weaker in 
the Alk-HN arrangement. To quanti fy th is we have per- 
formed MuUiken population analysis [SJI • The strongest 
effect is observed for LiNH, where the partial occupancy 
of the valence Pz orbital is 0.086 in the Li-NH arrange- 
ment compared to 0.004 in the Li-HN arrangement. The 
Pz occupancy for Alk-NH arrangements decreases down 
the periodic table: Na 0.068, K 0.037, Rb 0.027 and Cs 
0.021. This hybridization tendency is similar to that 
found for the alkali-metal trimcrs in the quartet states 

The present calculations give a quartet well depth for 
RbNH that is about 12% deeper than the value of 78 meV 
(630 cm~^) obtained by Soldan and Hutson using 
multireference configuration interaction (MRCI) calcula- 
tions. The present work used RCCSD(T) calculations, 
which give a better treatment of dispersion effects, and 
also used much larger basis sets. Our well depths are also 
substantially deeper than those reported for Rb-NH and 
Cs-NH by Tacconi et al. jSllusing methods and basis 
sets similar to those of Ref. [28[ . 

For all the Alk-NH systems there are also covalcnt 
states of symmetry (^A" at bent geometries) and 
an ion-pair state of ^H symmetry (^A' and ^A" at bent 
geometries). As described above, there is an avoided 
crossing between the two ^A" states. Table U gives the 
equilibrium positions and well depths of the ^H states 
for Alk-NH geometries, and it may be seen that for all 
systems the ion-pair well is more than 13000 cm~^ deep. 
The potential curve for the state cannot be obtained 



from RCCSD(T) calculations, but it is qualitatively simi- 
lar to that for the '^'E~ state in the long-range region [28| . 
Table [J includes the position and energy of the crossing 
point between the ^S" and ^H curves, and it may be seen 
that the crossing is always outside the minimum of the 
S state. Because of this, the lowest adiabatic surface of 
either ^A' or ^A" symmetry always has a very deep well 
of ion-pair character. This well is strongly anisotropic, 
so that any collision that samples the doublet surfaces is 
likely to be strongly inelastic. 

In conclusion, it appears that both the doublet 
and quartet states of Alk-NH systems have sufficient 
anisotropy to prevent sympathetic cooling for low-field- 
secking molecular states. 



B. Alkaline-earth atom + NH interactions 

The potential energy surfaces for alkaline-earth atoms 
(Ae) interacting with NH are substantially different from 
those for AlkNH systems. At linear geometries, there is 
one dispersion-bound state, which correlates with 

the Ae(^S') -I- NH(^I]") dissociation hmit. This state is 
crossed by ion-pair '^H and ^H states, which in a diabatic 
representation correlate with the Ae'*"(^S') + NH^(^H) 
dissociation limit. At non- linear geometries (point group 
Cs), the "^S" state becomes a ^A" state, and the "^H state 
is subject to the Renner- Teller effect and splits into two 
states with the electron hole either in the triatomic plane 
{^A') or perpendicular to it ('^A"). In cuts at fixed NH 
bond length, the dispersion-bound ^A" and ion-pair ^A" 
states form a conical intersection at linear geometries, 
while in the full three-dimensional picture they form a 
seam of conical intersections parameterized by the N-H 
distance. 

The counterpoise-corrected equilibrium distances and 
well depths for AcNH systems are shown in Table HI] for 
both Ac-NH and Ae-HN linear geometries. The corre- 
sponding potential curves arc shown in Figs. [2] and [3] for 
Ac-NH and Ae-HN geometries, respectively. It may be 
seen that the anisotropy of the dispersion-bound state is 
considerably smaller for the alkaline-earth atoms than for 
the alkali-metal atoms. However, for Ca and Sr the dif- 
ference between the well depths at the two linear geome- 
tries (60.4 and 184.6 cm~^, respectively) is still several 
times the NH rotational constant {B^ = 16 cm~^). For 
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Be and Mg, by contrast, the difference is only 10.5 cm"'^ 
and 3.6 cm~^, respectively, which is smaHer than the NH 
rotational constant. The difference may be understood 
in terms of the smaller s-p excitation energies for Ca and 
Sr (1.9 and 1.8 eV, respectively) compared to those of Be 
and Mg (both 2.7 eV). 

Another important feature is the position of the cross- 
ing between the dispersion-bound '^E" state and the ion- 
pair '^n states. These are also tabulated in Table HIl and 
shown in Figs.[2]and[3l For Be-NH and Mg-NH the cross- 
ing occurs fairly high on the repulsive wall of the "^S^ 
state, while for Ca-NH and Sr-NH it occurs at negative 
energies (in the potential well). This may be crucial for 
coUisional properties. If the crossing is located at nega- 
tive energies, as for Ca-NH and Sr-NH, the deep, strongly 
anisotropic ion-pair well may be accessed in low-energy 
collisions and is likely to result in strong inelasticity. On 
the other hand, if the crossing occurs at a high energy 
in a classically inaccessible region, as for Be-NH and Mg- 
NH, the deep ion-pair well is accessible only by tunneling 
through a barrier and may not have a strong effect on 
collisions. 

For Ae-HN geometries, shown in Table [Hi the cross- 
ing occurs high on the repulsive wall for all the AeNH 
systems. 



C. Mg-NH interaction potential 

As shown above, the BeNH and MgNH systems have 
potential energy surfaces that appear promising for sym- 
pathetic coo ling . However, Mg has been successfully 
laser-cooled [36[ whereas Be has not. We therefore fo- 
cus in this section on developing a complete potential 
energy surface for interaction of Mg with NH(^I]~). 

For bent geometries near the conical intersection, the 
two lowest triplet states of A" symmetry are near- 
degenerate. Under these circumstances single-reference 
coupled-cluster calculations are inappropriate. One al- 
ternative, which we have previously applied for RbOH 
[29I . [30I I , is to carry out multireference configuration in- 
teraction calculations including single and double exci- 
tations (MR-CISD). However, for MgNH the contribu- 
tion from triple excitations is extremely large: for exam- 
ple, the well depths of the linear state is underesti- 
mated by 40% in RCCSD calculations. Because of this, 
we use SAPT(DFT) calculations to study nonlinear con- 
figurations of MgNH. It has recently been demonstrated 
[50I . [U that SAPT(DFT) gives reasonably good resuhs 
for Mg2, NH-He and MgHe, and the polarizabilities and 
Van der Waals coefficients for Mg2 are reproduced with 
an accuracy of a few percent. A further advantage of the 
perturbation theory is that, by starting from zeroth-order 
wavefunctions corresponding to neutral monomers, we 
produce diabatic potential energy surfaces corresponding 
to neutral Mg-NH without contamination from ion-pair 
states. 

Fig. 13] shows the potential energy surface for Mg-NH 



obtained from SAPT(DFT) calculations. There are min- 
ima at both linear geometries and a saddle point between 
them, with a barrier of only 24 cm~^. The agreement 
between CCSD(T) and SAPT(DFT) methods for hnear 
configurations is very good for the Mg-NH geometry, 
where SAPT(DFT) overestimates the RCCSD(T) well 
depth by 6%. The agreement is slightly worse for the Mg- 
HN geometry, where SAPT(DFT) underestimates the 
weh depth by 9% (see Fig. [5]). 

The ion-pair state does not reach negative energies un- 
til distances R < 2.6 A. The dispersion-bound state is 
strongly repulsive at such distances. We have carried 
out multireference self-consistent field (MCSCF) calcula- 
tions of the two states of ^A" symmetry in the region of 
their avoided crossing for a range of angles using a cc- 
pVQZ basis set [4^|4y|. The lowest barrier for crossing 
onto the ion-pair state occurs for an angle ~ 110° at a 
distance of i? ~ 2.4 A and an energy of 4-2200 cm~^ with 
respect to the atom-molecule threshold. The singlet ion- 
pair state is about 700 cm^^ shallower than the triplet 
near its equilibrium geometry and will therefore cross the 
dispersion-bound state at even higher energies. 

The collision energies of importance to sympathetic 
cooling are in the range between 1 and 100 mK (0.06 
cm~^). At such energies we believe that the barrier sep- 
arating the wells of the ion-pair and dispersion-bound 
states is wide and high enough to neglect the conversion 
from MgNH to Mg+NH~ and to perform collision calcu- 
lations only on the dispersion-bound surface. 



IV. CONCLUSIONS 

In this paper we have presented an overview of the 
interaction potentials of alkali-metal and alkaline-earth 
atoms with NH molecules in their ground -^S" state. The 
interaction potentials of quartet states of AlkNH systems 
are strongly anisotropic, with deep wells at Alk-NH ge- 
ometries. The bonding in the well region involves strong 
mixing of the s and p orbitals of the alkali-metal atom 
and is thus covalcnt in nature. For geometries close to 
Alk-NH configurations the quartet states are crossed by 
ion-pair doublet states in the energetically accessible re- 
gion. Because of the presence of the ion-pair state, the 
lowest doublet adiabatic potential energy surface has a 
very deep well. The anisotropics for both doublet and 
quartet states are so large that it is unlikely that sympa- 
thetic cooling of NH by alkali-metal atoms will be suc- 
cessful for molecules in low-ficld-seeking states. 

For alkaline-earth atoms the interaction potentials are 
much shallower and less anisotropic, especially for BeNH 
and MgNH. The sp mixing is much weaker and the 
bonding is dominated by induction and dispersion forces. 
For MgNH the anisotropy is comparable to or smaller 
than the rotational constant of NH, B^- For CaNH and 
SrNH the anisotropy is rather larger, of the order of 
4Be for CaNH and IIB^ for SrNH. The relatively weak 
anisotropy raises the hope that alkaline earth atoms could 
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TABLE 11; Lowest and states of linear AeNH: minima {Rn 



ments (Ae-NH, Ae-HN). Energies are given in cm ^ and distances in A. 



and crossing points (_Rx, V^) at different arrange- 



Ae 






Ae-NH 










Ae-HN 


^min 


* mm 






nil 


^ mm 


-'•'min 


mm 


Be 


3.995 


-84.5 


2.30 


2390 


1.55 


-20240 


4.301 


-95.4 


Mg 


4.157 


-106.5" 


2.59 


1510 


1.95 


-10120 


4.636 


-103.0'' 


Ca 


3.963 


-165.7 


3.19 
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FIG. 2: One-dimensional cuts through the lowest (full) and '^H (dashed) potential energy surfaces of AeNH systems at the 
linear Ae-NH arrangement. For clarity, the region near the crossing of the dispersion-bound and ion-pair states is magnified. 



be used for sympathetic cooling of NH molecules. 

The dispersion-bound ^Y,^ states of AeNH systems are 
crossed by singlet and triplet ion-pair states and ^H. 
At nonlinear geometries the '^S" state becomes '^A" and 
there is a component of the •^11 state of the same sym- 
metry. The absolute minimum thus has ion-pair charac- 
ter in all cases. For CaNH and SrNH the ion-pair state 
crosses the dispersion-bound state in the energetically 
accessible region, so that the adiabatic potential energy 
surfaces have a single- minimum structure with a deep po- 
tential well and strong anisotropy. For BcNH and MgNH, 



however, the crossing occurs on the repulsive wall of the 
potential of the dispersion-bound state and the deep ion- 
pair well is likely to be inaccessible in low-energy colli- 
sions. We have calculated a full 2-dimensional potential 
energy surface for MgNH and verified that the crossing 
occurs on the repulsive wall at all geometries. 

The BeNH and MgNH systems are thus promising can- 
didates for sympathetic cooling. In future work we will 
carry out collision calculations on MgNH to explore this 
further. 
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FIG. 4: Contour plot of the potential energy surface for the 
dispersion-bound state of Mg-NH, from SAPT(DFT) calcula- 
tions (in cm~^). 
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